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INTRODUCTION 


Product random processes are used to model atmospheric turbulence in 
aeronautical applications. One example is the process developed by Press 
and his associates at NACA (references 1 and 2). This process is formed 
by the product of two Gaussian processes, one representing the local 
Gaussian properties of the turbulence and the other introducing a random 
amplitude modulation. This process is the basis of the current application 
of random process theory to the specification of aircraft strength 
requirements for flight in turbulence (references 3 and k) and has had 
extensive application to the analysis and measurement of aircraft response 
to atmospheric turbulence (references 5 through 7* for example). The 
mathematical properties of the Press model have been examined in reference 8. 
That study shows how the random process accounts for both the Gaussian 
property of short term turbulence experience and the non-Gaussian form of 
long term experience. In the present study the Press model is referred 
to as either a product or an amplitude modulated random process. 

Another random process which is used to model atmospheric turbulence 
is that developed recently by Reeves (references 9 through 12). This 
process is a modification of the product process by the addition of an 
independent mean value process, which allows more versatility in modelling 
the properties of measured turbulence data. The original development of 
references 11 and 12 was primarily concerned with the simulation of the 
short term properties of atmospheric turbulence. Consequently, the three 
components of the total process (called the local, amplitude, and mean 
value components in the present study) were specified to have the same 
integral scale value. 
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The present study is an extension and reinterpretation of the 
development of references 11 and 12. The basic random process is re- 
interpreted by introduction of the modulation concept: the local 

Gaussian component is modulated by the slowly varying amplitude component 
and by the addition of the slowly varying mean value component. The 
resulting random process is referred to as either the amplitude- 
mo dul at ed-plus-me an (AMPM) or the total process. The first part of the 
study is a general development of the first (at one time value), second 
(at two time values), and higher order properties of the process without 
restriction to the modulation condition. An approximate method, based on 
the modulation concept, for the analysis of the response of linear dynamic 
systems to the AMFM process is developed. Finally, the transition 
properties of the process, which explicitly show the modulation effects. 


axe examined. 



SYMBOLS 


A 

b 

C( ) 
c 

det( ) 

E [ ] 
erf ( ) 
erfc ( ) 
h (t) 
i 

M u 

m, M 
N ( ) 

N 

or 

P ( ) 

F (w | s) 

r, R 

s, S 
T 

t 

w, W 
z, Z 


standard deviation of either the subscripted process or 
(without subscript) the R process 
standard deviation of S process 
characteristic function of subscripted process 
standard deviation of M process 
determinant of matrix 
ensemble average 
error function (reference 13) 
complementary error function = 1 - erf ( ) 
impulsive response function of a linear system 
unit of imaginaries, J-l 
fourth order flatness factor (or kurtosis ) 
mean value process 

expected rate of positive slope crossings of indicated level 

expected rate of positive slope zero crossings of R 
1 A * 

process = — f 
r 

probability density function of subscripted process 
conditional probability density function of the v process 
conditional on the value of s 
local random process 
amplitude random process 
time interval 
time 

total or amplitude-modulat ed- plus-mean (AMPM) random process 
product or amplitude modulated random process 
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k 


a 


e 

p(t) 

a 

T 

$> (a>) 
4> (w) 

¥ (t) 
to 


ratio of the standard deviations of the Z and M 
processes 

Fourier transformation variable for characteristic function 
autocorrelation function of subscripted process 
modulus of amplitude random process 
time difference variable 

power spectral density function of subscripted process 
normalized power spectral density function of subscripted 
process 

autocovariance function of subscripted process 
frequency, Fourier transformation variable of T 


Superscripts : 


£ 


integer index 


Subscripts: 

c 

d 

i 

in 

i 

out 

a 


conditional 

derivative 

integer index 

input random process 

integer index 

output random process 

quasi-steady 


tr 


transition 


DEVELOPMENT OF THE AMPLITUDE-MODULATED-PLUS-MEAN (AMPM) RANDOM PROCESS 

The mathematical properties of the AMPM process are developed in the 
present section. The development starts with the definition of the AMPM 
process in terms of its three component processes. The probabilistic 
structure of the process , that is, the first order and the higher order 
properties, is then developed. No restriction is placed on the relative 
values of the integral scales of the three component processes. The 
development thus supplements that of references 9 through 12. 

Process Definition 

The product random process is defined by the product of two component 
processes. 


z(t) = r(t ) s(t) (1) 

The total or AMPM process is defined by the sum of the product process and 
a third component process. 


v(t) = z(t) + m(t) (2) 

v(t) = r(t ) s (t ) + m(t) (3) 

The three component processes are specified to be statistically independent. 
They are identified as: the local process R, the amplitude process S, 

and the mean value process M.* This terminology anticipates some 

•The present notation is based upon that of reference 8. The R, S, and M 
component processes correspond to the a(t), b(t), and c(t) processes, 
respectively, of references 9 snd 12. 


5 


6 


conditions which are subsequently placed on the relative values of the integral 
scales of the components: the rapidly varying local process R is modulated 

by a slower amplitude process S. The random process 2 will be referred 
to as either the product process or tbe amplitude modulated process. The 
random process W will be referred to as either the total process or the 
AMFM process. In both cases the latter terms will be used when the modulating 
properties of the amplitude component process are being emphasized. The 
former terms are used in the present general development, since the 
modulation condition is not introduced. 

First Order Properties 

The probabilistic structure of the total process is completely defined 
by the structures of the three component processes. The development starts 
with the formulation of the first order properties, that is, the random 
process is considered at one value of time only. The formulation is developed 
in terms of the characteristic functions, which are the Fourier transfor- 
mation of the probability density functions (reference lh). 

c v (e) = E[e iw0 3 (U) 

The characteristic function cf the product process is developed from its 
two components, equation (l), by using the relation for the product of two 
independent random processes (references 8 and ll+). 


c z (6) 85 C r (0s) p(s) ds 


(5) 
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The characteristic function of the total process follows from the defining 
relation, equation (2), which states that the process is the sum of two 
independent random processes. 

C w (0) * C z (0) C n (6) (6) 

Equations (5) and (6) express the characteristic function of the total 
process in terms of the characteristic and probability density functions 
of the three component processes. 

The probability density function of the total process similarly can be 
expressed in terms of the functions of the three component processes. 

The relationship is formulated by writing the Joint probability density of 
the total process and two of its components. This is developed by using 
the coordinate transformation of equation (3) and the independence of the 
component processes (reference 14). 

p(w,s,m) = -|~y P r P ( ® ) pU) (7) 

The probability density function of the total process is obtained from the 
joint function by integrating over the range of the two component processes, 
that is, as a marginal function of the joint density function. 

p(v) * P( w » s » m ) ds dm (8) 

Equations (7) and (8) are equivalent to equations (5) and (6). This 
equivalence can be shown by using tbe appropriate Fourier transformation 
relations (reference 15). 
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The total process can also he developed in terms of a conditional 
process 5 which is related to the local R component. The conditional 
process is the value of the total process for given values of the 
amplitude and mean value components. The associated conditional probability 
density function follows from the Joint density function „ equation (T ) 9 
and the definition of conditional probability, 

p c (w| Ss m) * ]^P r (— ^ (9) 

The preceding development has considered the characteristic and 
probability density functions of the component processes in general terms ; 
no specific form has been assumed for these processes. Following the 
developments of Press {references 1 and 2) and Reeves (references 9 through 
12 ) s the three independent component processes are specified to be stationary 
and Gaussian with zero mean values. The first order properties of the 
component processes are completely determined by their variances, which are 

E[r 2 ] = A 2 

E[s 2 ] = b 2 

E[m 2 ] * c 2 (10) 

This notation is guided by that used for the Press model of atmospheric 
turbulence in the aeronautical literature. * The probability density and 

*The notation of A 9 b a and c for the standard deviations of the component 
processes in the present study corresponds to that of a a9 c ^ 9 and C c 
respectively in reference 9 . 
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characteristic functions of the three component processes hare the Gaussian 
functional form. For the B component process the functions are 


' ■— i. 

. y 2 m A 

(11) 

- i e 2 A z 

c r (8) * e 

(12) 


The probability density and characteristic functions of the other 
components have similar functional foms; the notation for the associated 
variances is given in equation (10). < 

The statistical moments of the total process are developed from the 
defining relation, equation (3), and the moments of the three independent 
components. All odd order moments are zero. The relations for the variance 
and some of the higher order moments are 

E[v 2 ] = A \ 2 + c 2 

Etv* 4 ] = SA^b* 1 + 6A 2 b 2 c 2 + 3c 1 * 

E[v 6 ] = 225A 6 b 6 + 135 aVc 2 (13) 

+ U5 A tVc 1 * + 15c 6 

The relations for the moments can be expressed in an alternate notation. 
Since the total process is the sum of the independent product and mean value 
processes, the ratio of the standard deviations of these tvo processes i- a 
fundamental parameter of the total process. 
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a « Ab/c ' (1*0 

If the standard deviation of the product process approaches zero, then the 
value of a becomes zero. In this case the total process is dominated by 
the mean value process. If the standard deviation of the mean value 
process approaches zero, theu the value of a becomes infinite. In this 
case the total process is dominated by the product process. 

The statistical momeutr, can be expressed in terms of the a parameter 
by using equations (13) and (1*0 . 

E[v 2 ] = c 2 (a 2 + 1) 

E[v*] * 3c k (3<x k + 2a 2 + 1) (15) 

E[v 6 ] = 15c 6 (l5a 6 + 9a U + 3a 2 + 1) 


The associated fourth order flatness factor (or kurtosis) of the total 
process is 


Mj^(w) = 


Etv 1 *] = 3(3a^ + 2a 2 + l) 
E 2 [v 2 ] (a 2 + l ) 2 


(16) 


This factor shows the probabilistic structure of the total process in a 
concise manner. In the limit of small values of the a parameter, the 
total process becomes the Gaussian mean value process; the flatness factor 
has the Gaussian value of three. In the limit of large values of the a 
parameter, the total process becomes the product process, which is the 
product of two independent Gaussian processes. The resulting flatness 



n 


factor is equal to nine, which Is the square of the Gaussian value. 

The probability density and characteristic functions of the total 
process are obtained from the functions of the three component processes. 
The characteristic function of the product process is obtained by 
using equation (5) (references 8 and 11). 

c (e) ■ (aVb 2 + i )~ 1/2 (i?) 

Z 

The characteristic function of the total process is obtained by using 
equation (6). 

_ k . i e 2 c 2 

C (B) * (A^B 2 +1) 2 e 2 (l8) 

w 

The developsent of the corresponding probability density function involves 
an inverse Fourier transformation which appears to be intractable. The 
formulation can be extended to include the first order characteristic 
functions of a vector random process. However these functions are easily 
obtained from the general form of the higher order function, which is 
developed subsequently. 

In the original development (references 1 and 2) the product process 
was formulated in terms of the associated conditional process, equation 
(9)» rather than in terms of the local R component process. The 
probability density function of the conditional process associated with 
the total process is determined by equation (9) and by the specific form 
of the density function of the R component, equation (11). 
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P c (v|s,m) 


_ Iil-.bJI 

1 e 2A 2 s 2 

a| s | 


(19) 


Thus the conditional process is Gaussian with given values of the standard 
deviation and mean. In this approach the total process is formed from the 
conditional process by introducing the random variations of the standard 
deviation and mean value. In the present formulation the amplitude compo- 
nent is not the standard deviation of the conditional process, but can be 
related to it by giving the standard deviation factor A a uni t, value and by 
introducing a new random process which is restricted to nonnegative values. 


a = |sj 


( 20 ) 


The probability density function of the new process follows from that of 
the amplitude process by using the indicated transformation. Since only the 
modulus of the amplitude random variable appears in the probability density 
function of the conditional process, equation (19), and in the joint 
probability density function, equation (7), the formulation of the first 
order properties can be developed in terms of either of these two processes. 

Second and Higher Order Properties 

The formulation is extended to the higher order properties, that is, 
the relations for the total process considered at more than one value of 
time. The defining relation of the total process, equation (3), is valid 

v 

at all time alues. It defines the development in time of the total process 
in terms of the development of the three component processes. The higher 
order properties of the total process are thus defined by those of the 
three component processes. 



The autocovariance function is the joint second moment of the random 
process (with zero mean value) at two time values. Since the component 
processes are specified to be stationary, the autocovariance of the total 
process is a function of the difference between the two time values only. 

BW^) v(t 2 )] * ^(x) (21) 

where x = Jt g - t.J 

Using the defining relation, equation (3) , and the independence of the 
component processes, the autocovariance function of the total process is 
determined by the functions of the three components. 

yt) = y-o + yx) 

where (22) 

\(x) * yx) T s (t) 

The power spectral density function is the Fourier transformation of 
the autocovariance function. Using equation (22), the power spectral 
density function of the total process is expressed in terms of the 
functions of the component processes (reference 16 ). 

\(a) = $ 2 (a) + «ta(w) (23) 


where 
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The indicated integral is the convolution of the functions of the local and 
the amplitude component process. The convolution operation follovs from 
the Fourier transformation of the product of the two autocovariance functions 
in equation (22). 

The derivative of the total process is expressed in terms of the 
derivatives of the component processes by using the defining relation, 
equation (3). 

• • • • 

v(t) * r(t)s(t) + r(t)s(t) + m(t) (24) 

(This assumes the existence of the indicated derivatives.) The derivative 
is thus formed by the sum of two product processes and one mean value 
process. The autocovariance function of the derivative can be expressed 
in terms of the autocovariance functions of the components either from 
equation (24) or from the autocovariance relation, equation (22). 

The higher order probability density and characteristic functions of 
the total process are determined from the corresponding functions of the 
three component processes. The characteristic functions are generally 
easier to develop. The second order characteristic function of the 
product process is 

00 00 

C z^ 6 1* 6 2^ = ! f C r ( e :L S l* e 2 8 2^ P( s l» s 2^ dS l ds 2 (25) 

where s^ = s(t^) 

s 2 * s ^2^ 
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Since the total process is the sum of the independent product and mean 
value processes, equation (2), the characteristic function is the product 
of the functions of those tvo processes. 

C w (e i» 0 2 ) * C s (0 l* 0 2 ) C m (0 l* e 2 ) (26) 


The formulation can also be developed in terms of an associated 
conditional process. The relation for the second order probability 
density function of the conditional process is similar to that for the 
first order function, equation (9)* 


v. - 


1 ”1 " “1 w 2 “ m 2 
p c ( v v 2 lw b 2* b 2> “ TW P r ( ~ 


( 27 ) 


The associated second order characteristic function is obtained by taking 
the Fourier transformation of the probability density function. Using 
the appropriate identities or the Fourier transformation (reference 15 ), 
the relation is 


^ c ( 1 * ^ r (02 s i*®2 S 2^ 

ie^ + iQgffig 


( 28 ) 


Since the three independent components are specified to be Gaussian 
random processes, specific relations can be developed for the higher order 
characteristic functions of the total process. The properties of the 
Gaussian component processes are defined by their autocovariance functions, 
which are conveniently expressed in terms of the standard deviations and 
autocorrelat ion functions . 
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Ot) - aVp,(t) 

A A 

T (t) = c 2 P(t) (29) 

m m 

where 

P z (t) * P r (T)p g (x) 

(The tem autocorrelation function is used herein to identify the normalized 
autocovariance function. ) In a similar manner the power spectral density 
functions of the three component processes are defined in normalized form. 


$ (to) = A^b 2 ^ (w) 

z z 

p 

$ (to) = c 4> (to) 
m m 


( 30 ) 


The normalized power spectral density functions are the Fourier transfor- 
mations of the corresponding autocorrelation functions. The power spectral 
density functions are normalized in the sense that the integrals of the 
functions over the total frequency range have unit values. This property 
corresponds to that of the autocorrelation functions having a unit value 
when the time difference variable is set to zero. 

The autocorrelation and normalized power spectral density functions 
of the total process are related to the corresponding functions of the 
product and mean value processes. The relations are developed by 
combining equations (22) and (23) , which relate the autocovariance and 
spectral functions of the various random processes, with equations (29) 
and (30), which relate the total and normalized functions of the component 


processes 
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p v (T) “2 ^ 1 2 p m^ 

( 31 a) 

a + 1 a + 1 


2 . 


“2 , + 2 , <l> m (a)) 

(31b) 


These two relations show the contributions of the product and mean value 
processes to the functions of the total process. In the limit of small 
values of the a parameter, the functions of the total process become 
those of the mean value process. In the limit of large values of the a 
parameter, the functions of the total process become those of the 
product process. 

The expression for the second order characteristic function of the 
total process is developed from the characteristic and probability density 
functions of the three component processes. For example, the second order 
characteristic function of the E component process is 


c r (e 1# e 2 ) = exp{- ~a 2 [6^ + 2p r (x)e 1 e 2 + e|3> (32) 

The function is dependent upon the difference between the two values 
of time, which appears in the autocorrelation function. The second order 
characteristic functions of the other component processes have the same 
form, with the appropriate notation for the corresponding moments indicated 
in equation (29). 

The second order characteristic function of the total process is 
developed from equations (25) and (26) by using the appropriate functions 
for the components and performing the required integrations. 
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WV ■ 1(1 ' p r Kl - 0> 2 A * A* 2p r p s B l B 2 

+ b| + • exp[- | c 2 (e 2 + 2 p m e 1 e 2 + e|)] (33) 

where 

8i = Ab0 1 i = 1,2 

The associated second order . probability density function is obtained by 
the inverse Fourier transformation' of the characteristic function in both 
variables. The resulting inversions appear to be intractable. 

The second order characteristic function can be used to develop other 
quantities of interest. An example is the development of the 
characteristic function of the first derivative of the total process, 
which is obtained by using the approach of reference 17. 

C.(6 J = 11m C (- -4, -4) 
v d T -*> w T T 

- 1/2 - 1/2 

= [i - AVpjtoe®] [i - aVp^o)^] (3U) 

exp[“ c 2 p"(O)0^] 

The characteristic function of the derivative is the product of three 
separate characteristic functions. The first two of these are the 
characteristic functions of the product of two independent Gaussian processes, 
equation (17). The third 1? the characteristic function of a Gaussian 
process. Thus the derivati/e of the total process is the sum of three 
independent processes of the form indicated. This corresponds to the 
form of the derivative of the total process indicated by equation (2h). 
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The formulation can be extended to the development of higher order 
properties. The characteristic functions of higher order can be developed 
from the higher order functions of the Gaussian component processes. The 
following notation is used for the autocorrelation functions of the R 
component process. 

Elrft^rCtj)] * A 2 p r>iJ (35) 

Similar notation is used for the amplitude and mean value components. 

Using the independence of the product and the mean value processes, the 
characteristic function of general order is the product of the functions 
of these two processes. 

C;(5) = Cj(5)c-(S) ( 36 ) 

The characteristic function of the mean value component is that of a 
Gaussian process. The characteristic function of the product process is 
given in reference 8. 


where 


C~(6) = [det(q^ + 

q ij = kZb \ P r,i£ 8 fl. P s,£j 


(37 ) 


6 . , * Kronecker delta function 
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This relation is developed as the general order characteristic function 
of a scalar process, that is, for a scalar process at a general number 
(equal to the dimension n) of time points. The relation is also the 
joint characteristic function of a vector process at a single time point, 
with appropriate interpretation of the correlation coefficients of the 
component processes. Similarly the relation is the characteristic 
function of a vector process at several time points, with appropriate 
interpretation of the correlation coefficients. Thus equations (36) and 
(3T) are the general development of the characteristic functions of the 
total random process . 

QUASI-STEADY APPROXIMATION 

An approximate method for the analysis of the response of linear 
dynamic systems to the AMPM process is developed by using the modulation 
concept. This method, termed the quasi-steady approximation in reference 
8, was inherent in the original development of the amplitude modulated 
process (references 1 and 2). The resulting analytical technique has had 
extensive application in the analysis of aircraft dynamic response to 
atmospheric turbulence (references 5 and 6, for example). The concept of 
the quasi-steady approximation is applied to the AMPM process in the 
present section. 

The modulation concept introduces the idea of random processes with 
significantly different time scales. The amplitude component is assumed 
to be sufficiently slowly varying relative to both the local R process and 
the system dynamic properties that it does cot significantly influence 
the dynamic aspects of the system response. The system dynamic response 
to the amplitude modulated process is thus due to the rapidly varying R 
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component process; the amplitude component Influences the response only in 
a static manner. Since the dynamic response is entirely due to the local 
B process, vhich is Gaussian, the probabilistic structure of the response 
of the linear system can be completely determined. The addition of the 
mean value component introduces no fundamental difficulty in the analysis 
of linear system response, since that component is Gaussian and independent 
of the amplitude modulated process. 

The development of the quasi-steady approximation is shown by the 
analysis of the response of a linear dynamic system to the AMFM process. 

For linear dynamic systems the response or output is the convolution of 
the impulsive response of the system and the input function. 


w out (t) = I M-Ov^t - T)dT (38) 

o 

Using the defining relation, equation (3), the input process is expressed 
in terms of its three independent components. 


00 

v out (t) “ i h(T)r in (t ' T)s in (t * T)dT 


00 

+ / h(x)m. (t - T )dx 


(39) 


The system output process is formed by the two terms of equation (39)t 
vhich correspond to the system response to the amplitude modulated end 
the mean value terms of the input process, equation (2). The response 
of the linear system to the mean value component process can be 
analyzed exactly. 



Since the input mean value process is Gaussian, then so is the output 
mean value process. This is a general relation for linear systems; 
there is no restriction on the value of the integral scale of the mean 
value process. 

The quasi-steady approximation is applied to the analysis of the 
response of linear dynamic systems to the amplitude modulated process, 
which is the first term of equation (39). Two basic assumptions are 
made. First, it is assumed that the impulsive response function of the 
Bystem is essentially zero after some time period T, so that the function 
can be set to zero with negligible change in the computed system response. 
The upper limit on the time integral of the amplitude modulated term in 
equation (39) can accordingly be set to T. Second, the amplitude process 
is assumed to be sufficiently slowly varying, so that the amplitude is 
essentially constant over the system response period T. Since the 
amplitude is then constant during the integration period, it can be 
removed from the integration. Using these two assumptions, the system 
response, equation (39), becomes 


v . (t) s 
out 


■in (t) f l h(T)r in (t - T)dT 


m . (t ) 
out 




The integral, which now involves only the local R component process, is 
the same basic relation as equation (UO). 
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r out (t) ■ r l h(T> r in (t - T> dT (k2) 

The output R component process is Gaussian since it is the response of a 
linear system to a Gaussian process. 

The expression for the response of the linear system to the AMPM 
process is obtained by combining the previous relations. 

v out (t) * B in (t) + m out (t) (U3) 

Both the system input and output processes have the same probabilistic 
structure since they have the same defining relation, equations (3) and 
(43) .between the three independent and stationary Gaussian component 
processes. Thus the input and the system output are both AMPM random 
processes; they differ only in the forms of the autocovariance functions of 
the associated R and mean value component processes. The quasi-steady 
approximation is a method for developing the response of linear dynamic 
systems to the AMPM process. The restrictions of the method are that the 
local R component process is Gaussian and that the amplitude component 
process is much slower than both the R component process and all significant 
aspects of the system dynamic response. The mean value component process 
must be Gaussian*. In the present development all three component 
processes are stationary and Gaussian. 

•In the subsequent deve] nt the mean value component is also assumed 
to be slowly varying in We same sense as the amplitude component, 
although this property of the mean value component is not necessary 
for the validity of equation (43). 
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Using the quasi-steady approximation, the problem of the analysis of 
the response of linear dynamic systems is reduced to the determination of 
the variances (or the covariance functions) of the system response. The 
following notation is used for the variances of the components of the i. -put 
random process. 

* 4 ] - 1 

«[•“,] * t> 2 (MO 

* = 2 


The input R component process is defined to have unit variance. This 
condition establishes an arbitrary factor between the R and amplitude 
component processes, which is inherent in the defining product relation, 
equation (l). The following notation is used for the variances of the 
components of the output process. 


E[r 2 ] * A 2 
out r 


E[m out ] ■ & 2 


0*5) 


The relations for the variances of the system input and output process 
follow from the defining relation for the AMPM process, equation (3), 
and from equations (43) through (45). 
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E[v^J - b 2 + c 2 


Etv out ] » A ¥ + A a ' 2 


( 1 * 6 ) 


These relations express the variances of the two AMPM processes directly 
in terms of those of the component processes. Alternately, the variances 
can he expressed in terms of the variance of the mean value component 
process and the a parameter, which is the ratio of the standard deviations 
of the amplitude modulated and the mean value processes, equation (lk). 


E[vf n ] - c 2(4 ♦ 1) 


E[v 2 ] » A 2 c 2 (a 2 . + l) 
out m out 


( 1 * 7 ) 


The relation between the a parameters of the input and output 
processes is obtained by combining the previous relations. 


A b A 
_ r _ _r 

°out A c “ A a in 
m m 


( 1 * 8 ) 


This relation shows the effects of the system dynamics upon the 
probabilistic structure of the system output process. If the system is 
static, then the variances of the local H and mean value component 
processes are changed by the same factor; the a parameters for the 
input and the output processes are equal. The system input and output 
processes then have the same probability density functions, except for 
the difference in the values of the standard deviation of the total 
process. In the general case the dynamic properties of the system will 
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change the value of the a parameter. For a dynamic system which acts 
as a low-pass filter, the slower mean value component is increased relative 
to the R component process, resulting in a decreased value of the a 
parameter. For a system which acts as a high-pass filter, the R component 
is increased relative to the mean value component process, resulting in an 
increased value of the a parameter. Thus the system dynamic properties 
can change the relative contributions of the amplitude modulated and the 
mean value processes to the total random process, resulting in different 
functional forms of the probability density functions of the input and 
output processes. 

The relations for the quasi-steady >rm of the AMPM process are 
obtained from those for the exact process by setting the dynamic properties 
of the amplitude and mean value component processes to zero. The dynamic 
properties of the quasi-steady process are thus due to the R component 
process only. An example of this procedure is the development of the 
relations for the derivative of the quasi-steady process. The derivative 
of the exact process is the sum of three independent processes. 


w(t) - r(t) s(t) + r(t) s(t) 


+ m(t) 


(2k) 


The characteristic function of the derivative of the exact process, 
equation (3^), is accordingly the product of three characteristic functions. 
The characteristic function of the quasi-steady process is obtained from 
that for the exact process by setting derivatives of the autocorrelation 
functions of the amplitude and mean value component processes to zero. 
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c- (6 d } * El ' A 2 * 2 ^ 0 ) 6 ^ 1 (U9) 

Since this is the characteristic function for product of two independent 
Gaussian processes, the derivative of the quasi-steady process is 

v^t) = r(t) s(t) (50) 

The quasi-steady approximation has the effect of removing the mean value 
variation from the derivative of the AMPM process. The derivative 
operation acts as a high-pass filter which completely removes the mean 
value component, at least to the accvxaey of the quasi-steady’ approximation. 
For the quasi-steady derivative the ralue of the a parameter is 
^finite. The fourth order flatness factor of the derivative has a value 
of nine as opposed to the values between three and nine for the original 
process, as indicated by equation (l6). 

One important application of the quasi-steady approximation is the 
development of the associated exceedance expression, which is the 
expected frequency of positive slope crossings of a given level of the 
random process. The exceedance expression is developed from the joint 
probability density function of the random process and its first 
derivat ' (reference 18). Using the quasi-steady approximation, this 
joint probability density function is developed from that for two 
uncorrelated random variables at one time point. This in effect omits 
the local variation of the amplitude and mean value component processes. 

The resulting relation for the exceedance expression is 

00 00 

N(w) -If N(wjs,m) p(s) p(m) ds dm 

_«o 


( 51 ) 
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The exceedance expression of the conditional process is equal to the product 
of one-half of the first absolute moment of the derivative and the 
Gaussian conditional probability density function (references 18 and 19) • 

B(v|s^) = -^p (52) 

Y2rf 


Using the probability density functions of the three component processes 
and evaluating the integrals in equation ( 51 ) , the exceedance expression 
of the quasi-steady process is 


N(v) _ 1 e c 2 /2A 2 b 2 
K ”2 


or 


v 

“ Ab 

[1 + 

v2c rZAb 


v 

+ e Ab erfc(-~- + —— )) (53) 

Y2c *2Ab 


where 


w > 0 


H(-v) = N(w) 


This is expressed in terns of the standard deviations of the three component 
processes. An alternate form for the exceedance expression is 

1 

"or 2 <£c /Sa 

V 

+ e ac erfc(— + — )} (5*0 

fic S2a 
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vhere 

v > 0 

N(-v) *» N(v) 

This is expressed in terms of the standard deviation of the mean value 
component and the a parameter* The exceedance expression contains 
both the exponential dependence of the amplitude modulated process and 
the Gaussian dependence, through the error functions, of the mean value 
component process. The general case is a combination of these tvo 
functional forms. The expected number of zero crossings of the quasi- 
steady process is 

1 

2 

H(v * 0) * N e 2a erfc(-^-) < N (55) 

or fia ~ or 

It is thus important to distinguish between the expected number of zero 
crossings of the AMFM process and those of the R component process, 
since these tvo quantities are generally not equal. The difference 
between: the tvo is due to the quasi-steady effect of the mean value 
component . 

The expected number of zero crossings, equation (55), has tvo 
limiting cases. In the limit of a vanishing contribution from the 
mean value component, a becomes infinite; the expected number of zero 
crossings is equal those of the R component process. In the limit of a 
vanishing contribution from the amplitude modulated process, O' becomes 
zero; the expected number of zero crossings approaches zero. 

lim N(w = 0) = N or 
a*o 


f • 


( 56 ) 
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In this case the mean value component becomes dominant. However the 
expected number of zero crossings of the AMPM process does not approach 
that of the mean value process, since the dynamic properties of that 
component have been eliminated by the quasi-steady approximation in the 
development of the exceedance expression, equations (53) and (5^). 

Several limiting cases of the quasi-steady exceedance expression are 
of interest. In the limiting case of vanishing contribution from the mean 
value component, a becomes infinite; the exceedance expression reduces 
to the exponential form of the amplitude modulated process. 

lim N(v) = N e "l W ^ Ab (57) 

or 
<r*« 


This expression was the basis of the original development of the amplitude 
modulated process in reference 1. The limit of vanishing contribution 
from the amplitude modulated process (a becomes zero) is more complicated 
since the dynamics of the dominant mean value component have been eliminated 
by the quasi-steady approximation. However the general case of non-zero a 
and finite w can be considered. For small values of both the a parameter 
and the level of the AMPM random process, the exceedence expression 


approaches the Gaussian form of the mean value component process. 



N( |v|«c; a«l) 



( 56 ) 


For extreme values of the AMPM process , the exceedance expression approaches 
the exponential form of the amplitude modulated process regardless of the 
(non-zero) value of the a parameter. 
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1 

N(|w| » ca) « N e 2a e “l w |/ ctc (59) 

Thus if there is any contribution from the amplitude modulated process to 
the AMPM process, that is, if the Qt parameter is not exactly zero, 
the contribution from the amplitude modulated process vill dominate the 
exceedance expression at sufficiently large values of the AMPM random 
process. 

The preceding development of the exceedance expression is based upon 
the quasi-steady approximation, which eliminates the dynamic properties 
of the amplitude and mean value component processes. The validity of the 
resulting qua si- steady exceedance expression is questionable in the case 
of small values of both the a parameter and the level of the AMPM process. 
In this case it is necessary to develop the exact form of the exceedance 
expression of the AMPM process, following the approach of references 11 
and 12. The resulting exceedance expression apparently can not be 
expressed in analytical fora, but requires numerical integration procedures. 

TRANSITION PROPERTIES 

The transition properties associated with a random process show the 
nonstationery development of the process from an initial value of time. 

In the case of the AMPM process, the transition properties explicitly show 
the effects of the slowly varying amplitude and mean value components upon 
the total process. The transition properties show the development of the 
random process from the initial Gaussian fora of the local R component 
process with a given amplitude and mean value to the non-Gaussian form of 
the fully developed AMPM process. 
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The primary quantity of interest is the transition probability density 
function of the AMPM process at one time value, conditional on the value 
of the process at an earlier time value (references 19 and 20). For the 
AMPM process the transition probability density function is conditional 
on the earlier value of all three component processes. 

p[w(t 2 )|r(t 1 ), s(t x ), m(t x )] 

= ptvgjr^ s 1# m 1 3 t x < t 2 ( 60 ) 

The transition of the AMPM process between the two values of time occurs 
in two stages, since the component processes are classified as rapidly 
varying (the R component) and as slowly varying (the S and M components). 

In the first transition stage the density function becomes independent of 
the initial value of the local R process. This transition is not important 
in most applications and is emitted in the subsequent discussion. In the 
second transition stage the probability density function becomes independent 
of the initial values of the amplitude and mean value processes, thus 
explicitly showing the effects of the development of these two components. 

The transition probability density function is obtained from the 
conditional function of the original process and the transition functions 
of the amplitude and mean value components. 

00 00 

p( v 2 |s 1 ,mi) * / / p(w 2 , s 2 , n 2 |s 1 , n^) ds 2 dm g (6la) 

—00 —00 

00 00 

p( v 2 l s i»®x^ ■ / / p(v 2 js 1 , EL ; s 2 , nu) pfsgjs,) (6lb) 

—00 —00 


p(m 2 |a 1 ) ds 2 dmg 
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The first relation is the development of the marginal function from the 
joint transition probability density function. The second relation is 
the combination of the first vith the definition of the conditional 
probability density function. 

Equation (6lb) is simplified by using a functional property of the 
conditional probability density of the original process. 

p(v 2 |b 1 , m^; s 2 , m 2 ) = p(w 2 |s 2 , m^) (62) 

Thus the conditional form of the AMPM process, that is, conditional on the 
values of the amplitude and mean value components at the same time value, 
is independent of the values of those two components at an earlier time 
value. This relation follows from the functional properties of the 
associated characteristic function. 


C^gls^ m^Sg, m 2 ) * C(0 1 ,e 2 |s, Q (63a) 


i6 nn + i0 5 m 

= C r^ 0 l S l* 6 2 B 2^ e le x = 0 


(63b) 


The first equation expresses the characteristic function as a marginal 
function of the second order characteristic function. The second equation 
is the combination of the first with the functional relation, equation (28), 
for the Joint conditional characteristic function of the AMPM process. 

This equation shows that the conditional characteristic function and thus 
the conditional probability density function is independent of s 1 and 
m^, thereby establishing equation (62). 
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The relation for the transition probability density function is 
obtained by combining equations (6lb) and (62). 

00 00 

P( v 2 l®l*»l) * / / p(v 2 js 2 ,m 2 ) p(s 2 |s 1 ) 

■.OO mOO 

p(m 2 |m 1 ) ds 2 dm 2 (6h) 


The transition probability density thus can be developed from the first 
order conditional density function of the original process and the 
transition functions of the amplitude and mean value components. These 
functions are known from the basic formulation of the AMPM process. Since 
the component processes are stationary and Gaussian, their transition 
distributions are known from the associated mean and variance. For the 
amplitude component for example, the transition probability density 
function is 


p(s 2 |s 1 ) = 


szaj i - p 2 (x) b 


( 65 ) 


-1*2 - P„(T) s ] 

exp { — p > 

2b 2 [l - p?(x)] 


Since the time dependence of the transition functions of the amplitude 
and mean value components appears only in the associated autocorrelation 
functions, the development of the transition probability density function 
of the AMFM process, equation (61*), is dependent solely upon the 
autocorrelation functions of those two components. 
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The previous relations can he used to develop the limiting cases of 
the transition probability density function. The limiting functions of 
the amplitude component process are 

lim pUgls^ * 6(sg - s L ) 

lim p(s 2 |s 1 ) = p( Sg ) 

'f' >00 

where 6(s) = Dirac delta function 
Thus in the limit of small time differences the amplitude component is 
equal to its initial value s^ with probability one. In the limit of 
large time differences the amplitude component process reaches its fully 
developed form with the probability density function being independent 
of the initial value. Similar relations hold for the mean value component. 
The corresponding limiting cases of the transition probability density 
function of the total process are developed from equations (6M, (66a), 
and (66b). 


(66a) 

(66b) 


l v i “ “l 

liM p(» 2 |. 1 ,« 1 ) = » ]^JP r {— ) 


(67a) 


lim p(v 2 is ,m x ) = p(v ) (6Tb) 

X"W 


Equation (9) has also been used in the first relaticn. In the limit of 
small time differences, the transition function has a Gaussian form with 
given emplitude and mean value. In the limit of large time values, the 
transition function approaches the form of the fully developed process, 
which is independent of the initial values of the amplitude and mean value 


components. 
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The transition properties of the AMPM process can also he developed 
in terms of an associated transition random process. This nonstationary 
process corresponds to the transition probability density function, 
equation (6^). The transition process is formed from the fully developed 
R process and the transitional forms of the amplitude and mean value 
components. 

V r S tr + «tr ' t68) 

The definition of the transition process in terms of the associated three 
compoi ent processes follows from both equation (6lb) and the form of the 
relations for the transition process, which correspond to the appropriate 
relations for the original AMPM process, equations (7) through (9). 

The transition processes of the amplitude and mean value components are 
defined by their transition probability density functions, for example, 
equation (65). 

The transition process has an associated set of transition moments. 

00 

* / Wg p(v 2 js 1 ,m 1 ) dv 2 (69) 

_00 

The transition moments are obtained from the definition of the transition 
process, equation (68), by expressing the transition moments in terms of 
the moments of the independent component processes. The first four 
transition moments are listed in table 1. The transition moments are 
functions of time since they are functions of the autocorrelation 
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functions of the amplitude and mean value components* The transition 
properties allow the possibility of a non-zero values of the odd order 
moment s , which are zero for the fully developed form of the AMPM process. 

The transition moments show the two limiting cases of the AMPM proses® for 
small and large time values. In the limit of small time values, the mean 
value of the transition process approaches the initial value of the mean 
value component. In the limit of large time values, the transition moments 
approach the moments of the fully developed process, which are independent 
of the initial values of the amplitude and mean value components. All 
odd order transition moments become zero in the limit of large time values. 

The transition properties can also be examined by using the 
associated flatness factors, which are obtained from the transition moments. 
In the limit of small time values, the fourth order flatness factor is 


lim 

T+o 



/ .2 2 
(A s 




(70) 


This factor has a minimum value of one and a maximum value of three. In 
the former case the moments are dominated by the initial value of the mean 
value component. In the latter case the initial value of the amplitude 
modulation process is dominant; the flatness factor has the Gaussian value 
of three, which results from the fully developed R component process. In 
the limit of large time values, the fourth order flatness factor becomes 
that of the fully developed AMPM process, equation (16). 
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TABLE 1. Transition Moments of the AMPM Process 
E Kr ] ‘ 

« <,] * U - »y ♦ * 2 A 

E[, tr 1 * Vl 
E£*| r 3 * (l •• p 2 )* 2 + 

E[w tr ] - Efc tr ] 

E[v 2 r ) - A £ E[ S 2 r ] + E[m 2 r ] 

Etv^] « {3A 2 E(s 2 r l ♦ 3 E[m 2 J 

-2E 2 [m tr U E[m tr l 

E[vJ r ] » SaVIs^J-SaVIb^] + 6A 2 E[s 2 r ] E[m 2 r ] 

+ 3E 2 [m| r ] - 2E* , [in tr ] 


39 


The transition random process has an associated exceedance expression, 
vhich also shows the nonstationary effects of the development of the 
amplitude and mean value components. The transition exceedance expression 
is developed from the relations for the second order transition probability 
density function, which is similar to equation (6k) t from which the joint 
density function of the transition process and its first derivative can be 
developed. Using this Joint probability density function and the quasi- 
steady approximation, the relation for the transition exceedance expression is 

00 

N(w 2 |s 1 ,m 1 ) « If N(vJs 2 ,m 2 ) p(s 2 !s,) 

—00 

( 71 ) 

pOngjn^) ds 2 dm g 

All of the required quantities are known: the conditional exceedance 

expression follows from equation (52) and from the conditional density 
functions of the amplitude and mean value processes. Thus the transition 
exceedance expression can be determined. However it can not be obtained 
in explicit form since the double integration of equation (71) is apparently 
intractable. The transition exceedance expression shovs the limiting 
cases of the transition process. In the limit of small time values, the 
exceedance expression becomes that of the quasi-steady conditional process 
at a single time value, equation (52). In the limit of large time values, 
the exceedance expression becomes that of the fully developed process, 
equations (53) and (5*0* 


SUMMARY OF RESULTS 


The mathematical properties of a random process formed by the product 
of a local and an amplitude process, and the sum of that product with a 
mean value process are examined. The probabilistic structure of the re- 
sulting amplitude-modulated-plus-mean (AMPM) process is developed from the 
structures of the three component processes (the local, amplitude, and 
mean value components), which are specified to be stationary and Gaussian. 
The development includes the first and second order properties, and the 
characteristic function of general order. 

The quasi-steady form of the AMPM process is derived by omitting the 
dynamic properties of the amplitude and mean value components. A method 
for the analysis of the response of linear dynamic systems to the quasi- 
steady form of the AMPM process is developed. The system response is 
shown to be an AMPM process also, at least to the accuracy of the quasi- 
steady approximation. The probabilistic structure of the system response 
can however be different from that of the excitation process, depending 
on the properties of the dynemic system. An analytical relation is 
developed for the exceedance expression of the quasi-steady form of the 
AKPK process. The exceedance expression shows a combination of the 
exponential form of the amplitude modulated process and the Gaussian fora 
of the mean value component process. 

The transition properties of the AMFK process are examined. The 
transition properties show the nonstationary aspects of the process, which 
are due to the development in time of the amplitude and mean value component 
processes. The basic relations for the transition probability density and 
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characteristic functions are derived. The transition properties shov the 
development of the AMPM random process from the initial Gaussian form of 
the process vith given values of the amplitude and mean value to the 
non-Gaussian form of the fully developed process. 
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